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Abstract. Motivated by the problem of dealing with incomplete or imprecise 
acquisition of data in computer vision and computer graphics, we extend re- 
sults concerning the stability of persistent homology with respect to function 
perturbations to results concerning the stability with respect to domain per- 
turbations. Domain perturbations can be measured in a number of different 
ways. An important method to compare domains is the HausdorfE distance. We 
show that by encoding sets using the distance function, the multidimensional 
matching distance between rank invariants of persistent homology groups is 
always upperly bounded by the Hausdorff distance between sets. Moreover, 
we prove that our construction maintains information about the original set. 

Other well known methods to compare sets are considered, such as the sym- 
metric difference distance between classical sets and the sup-distance between 
fuzzy sets. Also in these cases we present results stating that the multidi- 
mensional matching distance between rank invariants of persistent homology 
groups is upperly bounded by these distances. 

An experiment showing the potential of our approach concludes the paper. 



Introduction 

Persistent topology is a theory for studying objects related to computer vision 
and computer graphics, by analyzing the qualitative and quantitative behavior of 
real- valued functions defined over topological spaces and measuring the shape prop- 
erties of the topological space under study (e.g., roundness, elongation, bumpiness, 
color) . More precisely, persistent topology studies the sequence of nested lower level 
sets of the considered measuring functions and encodes at which scale a topological 
feature (e.g., a connected component, a tunnel, a void) is created, and when it 
is annihilated along this filtration. At the very beginning of the development of 
persistent topology, this encoding captured only the connected component changes 
in the lower level sets of a real valued function, and took the name of size function 
[13 1 Some years later, it was extended to consider all homotopy groups of the 
lower level sets of a vector- valued function, under the name of size homotopy groups 
[14) . Nowadays we have a wide choice of variants for this encoding, ranging from 
persistent homology groups capturing the homology of a one-parameter increasing 
family of spaces [12] , to multidimensional persistent homology groups extending the 
previous concept to a multi-parameter setting [i] , to vineyards coping with changes 
in the function over time [7], to interval persistence [5], just to cite a few. In this 
paper we focus on multidimensional persistent homology groups. For application 
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purposes, these groups are further encoded by considering only their rank, yielding 
to a parametrized version of Betti numbers, called rank invariants (or persistent 
Bctti numbers). 

The stability of multidimensional rank invariants is quite an important issue 
in persistent homology theory and its applications because the lack of stability 
would make this invariant useless, every data measurement being affected by noise. 
Stability with respect to perturbations of the measuring function was proved in 
[5], based on the results of [U [5], comparing persistent homology groups by the 
multidimensional matching distance. 

In this paper we consider the problem of stability with respect to changes of the 
topological space, which is as much important as the stability with respect to the 
change of measuring functions. Changes of the space under study can be measured 
in a number of different ways. Indeed, according to the kind of noise producing 
the perturbation, some distances are more suitable than other to compare sets. For 
example, the Hausdorff distance is useful to measure distortions of the domain, 
while the symmetric difference distance can cope with the presence of outliers. 
Due to the existence of many different ways to compare sets, we propose a general 
approach to the problem of stability of persistent homology groups with respect 
domain perturbations, and we apply this approach in a few cases. 

Our main idea is to reduce the problem of stability with respect to changes of 
the topological space to that of stability with respect to changes of the measuring 
functions. This is achieved by substituting the domain K we are interested in 
with an appropriate function Jk defined on a fixed set D containing A', so that 
the perturbation of the set K becomes a perturbation of the function fx- As a 
consequence, the original measuring function ■ K — > M*' is replaced by a new 
measuring function $ : ^ K'^+i, $ = (/k, ip). Rank invariants of (D, $) can be 
compared using the multidimensional matching distance. In this way we can prove 
robustness of persistent homology groups under domain perturbations. 

In particular, we use this strategy when sets are compared by the Hausdorff 
distance. In this case, taking fx equal to the distance function from K , we prove 
that the multidimensional matching distance between the rank invariants associated 
with two compact sets Ki and K2 is always upperly bounded by the Hausdorff^ 
distance between Ki and K2 (Theorem 12. At the same time, we show that, 
in our approach, the information about the original domain K and its original 
measuring function if is fully maintained in the persistent homology groups of 
{D,S) (Theorem [13). 

As a further contribution, we show stability with respect to perturbations of 
the domain measured using the symmetric difference distance. Also in this case, 
associating with K a suitable function fx, we can prove stability with respect to 
domain perturbations measured by the symmetric difference distance (Theorem 
[3T|). 

We also consider the situation where sets are described in a fuzzy sense, by means 
of probability density functions, easily obtaining a stability result also in this case 
(Theorem [321) ■ 

An experiment on a binary image concludes the paper, illustrating our results. 

To conclude this introduction, we emphasize three key-points in this paper. 
Firstly, there is a real need for developing techniques providing persistent homology 
with robustness against domain perturbations since the classical setting is not stable 
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in this respect, as the example in Section |4] clearly shows. Secondly, the technique 
we propose in order to achieve stability makes sense only in the multidimensional 
version of persistent homology, since it is based on passing from a fc-dimensional 
measuring function to a (fc + l)-dimensional one. Thirdly, we underline that, sub- 
stituting the study of the domain K with that of a much simpler domain D, we 
obtain the searched for stability but not at the price of forgetting the persistent 
homology of K. 

1. Preliminaries 

1.1. Multidimensional persistent homology groups. The following relations 
^ and -< arc defined in R'"': for u = (iti,...,Ufc) and v = {vi, . . . ,Vk), we say 
u :< V (resp. it ^ v) if and only if Ui < Vi (resp. Ui < Vi) for every index i = 
1, . . . ,k. Moreover, M*^ is endowed with the usual max- norm: || (wi, 1*2, ■ • ■ , ^*fc)||oo = 
maxi<i<fc \ui\. 

We shall use the following notations: A+ will be the open set {(m, v) E M'"' x R'^ : 
u -< v}. Given a topological space X, for every fc-tuple u = {ui, . . . ,Uk) S MJ' 
and for every continuous function : X ^ R'^, we shaU denote by X{(p ^ u) the 
lower level set {x £ X : ipi{x) < Ui, i = 1, . . . ,k} and by ||<^||oo the sup-norm of 
if, i.e. Il'/'lloo = maxajgx ||(/?(a;)||oo. The function ip will be called a fc-dimensional 
measuring (or filtering) function. 

Definition 1.1. Let tt^"'^^ : IIq{X{(f d: u)) ^ Hq{X{0 < v)) he the homomor- 
phism induced by the inclusion map tt^"''") : X{0 < u) ^ X{ip < v) with u ^ 1/, 
where Hq denotes the gth Cech homology group. If u -< u, the image of tt^"'^'' is 
called the multidimensional qth persistent homology group of {X,ip) at (u, u), and 
is denoted by Hq^'"\x, ip). 

In other words, the group Hq^'^'^ {X, 0) contains all and only the homology classes 
of g-cycles born before u and still alive at v. 

In what follows, we shall work with coefficients in a field K, so that homology 
groups are vector spaces, and hence torsion-free. Therefore, they can be completely 
described by their rank, leading to the following definition (cf. [4]). 

Definition 1.2 {qth rank invariant). Let X be a topological space, and : X ^M.^ 
a continuous function. Let q £ 1. The qth rank invariant of the pair (X, 0) is the 
function P[x.0),q ■ — > N U {00} defined as 

P(x,0),q{u,v) = rankimTrf"'""). 

If X is a triangulable space embedded in some R", then P{x,ip),qiu, v) < +00, for 
every (u, v) e A+ and every q G Z [31 [5] . 

1.2. Matching distance. We now recall the construction of the distance Dmatch 
to compare rank invariants of multidimensional persistent homology groups. 

Definition 1.3 (Admissible pairs). For every vector I — {li,...,lk) of R'^ such 
that > for z = 1, . . . , fc and X)i=i — Ij ^.nd for every vector b= (61, . . . , 6fc) of 
R*^ such that Yli=i — 0; shall say that the pair (/, h) is admissible. We shall 
denote the set of all admissible pairs in R'^ x R'^ by Admk- Given an admissible 
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pair {l,b), we define the half-plane '''(f of MJ^ x M'^ by the following parametric 
equations: 




for s, t e M, with s <t. 

The key property of this foliation of A"*" by half-planes is that the restriction of 
P{x.0),q to each leaf can be seen as the 1-dimensional rank invariant associated with 
a suitable pair {X,F^^^^), where P^^f^j : X ^ M.. Precisely, the following statement 
(proved in holds: 

Theorem 1.4 (Reduction Theorem). Let {I, b) be an admissible pair and let F^-^ : 
X — > M &e defined by setting 

t,-,^ (x) = max < > . 

(l.b)"- ' ,= l,...,fc\ J 

Then, for every (u, v) — {si + b,tl + b) G '''(f J) we have that 

P(X,<p),q{^,v) = P(X F-^ ) ,(*'^) • 

Since 1-dimensional rank invariants can be compared by a distance dmatch that 
matches points of persistence diagrams [SJ |5] , we can construct the following dis- 
tance between fc-dimcnsional rank invariants, called multidimensional matching dis- 
tance: 

Drnatch{.P{X,(p),q, Pix tf) q) ^ ' '^match I P,x ) g' Pix G'* ) o ) ' 

(T.b)£Admk ' \ ' ' ' ii^.iy / 

where F^"^-^ = maxi=i^...,fc j '^'t^)-''' | and Gj^-- = max,=i,...,fe |M£h:^|, 

The key property of Dmatch is its stability with respect to perturbations of the 
measuring function ip, as the following theorem states (cf. [S]). 

Theorem 1.5 (Multidimensional Stability Theorem). Let X be triangulable. For 
every q E 1^, and for every two continuous functions (p,ip : X M'^ , 

Dmatch {p{X,0),q,Pt^x,-$).^ - \W ^ M\oo- 

1.3. Comparison of sets. The problems of description and comparison of sets 
can been dealt with in a myriad of different ways, each one more or less suitable 
than another for a given application task. 

In classical set theory, the membership of elements in a set is assessed in binary 
terms according to a bivalent condition - an element either belongs or does not 
belong to the set. By contrast, in fuzzy set theory [T71 [TT], a fuzzy set A m X 
is characterized by a membership function fA'-X^ [0, 1], with the value fA{x) 
representing the grade of membership of a; in A. Usually, the nearer the value of 
fAix) to 1, the higher the grade of membership of x in A. The fuzzy set theory can 
be used in a wide range of domains in which information is incomplete or imprecise. 

If classical set theory is adopted, then a number of different dissimilarity mea- 
sures exist to compare two sets [13 [TO] • A frequently used dissimilarity measure is 
the HausdorfJ distance, which is defined for arbitrary non-empty compact subsets 
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Ki,K2 of W\ If Ki,K2 are contained in a non-empty compact subset D of M", 
the Hausdorff distance can be defined by 

Sh{Kii K2) = maxjmax (a;), max cIkz [y]}, 

where dK denotes the distance to K, that is the function dx : D ^ M. defined by 
dif (a;) = min.ygif ||a: — y\\, \\ ■ j| being any norm on M" (c-g-, the Eucfidean norm). 
This can be reformulated as follows (cf. [SJ Ch.4, Sect. 2.2]): 

(1.1) 6H{Ki,K2)^\\dK,~dK,\U 

The Hausdorff distance is robust against small deformations, but it is sensitive to 
outliers: a single far-away noise point drastically increases the Hausdorff distance. 

A dissimilarity measure that is based on the area of the symmetric difference, 
such as the symmetric difference pseudo-metric, overcomes the problem of outliers. 
Denoting by n the Lebesgue measure on M", the symmetric difference pseudo-metric 
is defined between two measurable sets A^B with finite measure by d/\{A,B) = 
H{AAB) where AAB ^ {AU B)\{An B) is the symmetric difference of A and B. 
It holds that dA{A,B) = if and only if A and B are equal almost everywhere. 
Identifying two sets A and B if iJ,{AAB) = 0, we obtain the symmetric difference 
metric. 

Other dissimilarity measures for more restricted patterns are, for example, the 
bottleneck distance between finite point sets and the Frechet distance between 
curves. However, since many other distances could be considered, we will limit our 
research to consider stability with respect to the Hausdorff and symmetric difference 
distances. 

When fuzzy sets are used, their dissimilarity can be measured by any function 
distance. In this case we will confine ourselves to consider the sup-norm between 
fuzzy sets. 

2. Stability with respect to Hausdorff distance 

Our main idea in proving stability of rank invariants with respect to noisy do- 
mains is to transform perturbations of sets into perturbations of functions. In this 
way it is possible to apply the Multidimensional Stability Theorem 11.51 When do- 
main perturbations are measured by the Hausdorff distance, in order to pass from 
a set to a function, we insert the distance function da described in subsection 
11.31 as the first component of the measuring function. In this way, assuming that 
all the sets under study are contained in a larger set Z?, the original problem, i.e. 
studying persistent homology groups of a set K endowed with the restriction to 
K oi a measuring function ip : D ^ M.'^ , is transformed into the new problem of 
studying the persistent homology groups of D endowed with the measuring function 
$ = idK,<f). 

Given two domains Ki and K2, and two functions ipi,(p2 ■ D ^ M'^, our first 
result relates the distance Dmatch between the new pairs (D,$i), (Z?, $2) to the 
change of the measuring functions 01 and if 2, and to the Hausdorff distance between 
the original sets Ki, K2- More precisely, it proves stability with respect to both set 
and function perturbations. Indeed, the change in the multidimensional matching 
distance Dmatch is shown to be never greater than the maximum among the change 
in the Hausdorff distance between the domains Ki and K2 and the change in the 
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sup-norm between the measuring functions ipi and 02- In particular, if ipi and 
ip2 coincide then the change in the multidimensional matching distance Dmatch is 
never greater than the Hausdorff distance between Ki and K2. 

Theorem 2.1. Let Ki,K2 be non-empty closed subsets of a triangulable subspace 
D of M" . Let , rfif2 : -D — > K &e their respective distance functions. Moreover, 
let 01^02 ■ D R*^ be vector-valued continuous functions. Then, defining $1, $2 ■ 
D — > M.^^^ by "I>i = (dif 1 , ^i) and $2 — {dx^i 02), ths following inequality holds: 

Dmatch < Uiax (i^i , Xj) , \\01 " !^2||oo} ■ 

Proof. The Multidimensional Stability Theorem II .51 for measuring function pertur- 
bations implies that Dmatch {p[D,ii).q^ P[D,i2),q) - ^ *J'2||oo- It follows that 

Dmatch [P{D,$i),q^P{D,$2),q) - max{||dA'i - dK^Woo, \\01 - <^2||oo}- 

Hence, by equality (1.1), the claim is proved. □ 

We now consider the problem of retrieving the rank invariants of (X, 0\k) from 
the rank invariants of {D, $), with <f> = (dx, 0). The next result shows that for any 
sufficiently small value of /3 S M there exists a sufficiently small value a G R with 
< a < ^ such that P(if,<^|^),g(w, t?) = P(D,$),q W,v)). 

Theorem 2.2. Let K be a non-empty triangulable subset of a triangulable subspace 
D of R" . Moreover, let : D ^ R'^ be a continuous function. Setting $ : 13 — > 
jjfc+i^ (j) = {dx, 0), for every u,v £ R'^ with u < v, there exists a real number /3 > 
such that, for any /3 € M with < /3 < /?, there exists a real number a = Q!(/3), with 
< a < /?, for which 

P(K,0^k)A'^^ ^) = P{D,$),q ^)) ' 

for every a G R with Q < a < a. Ln particular, 

PiK,0^^),qiu,v) = lim^p^^^^^^^ ((0,u), (/3,iT)) . 

Proof. For every u <E R*^ , we have 

K{0\K di u) = {x Cz K : 0{x) ^ u\ 

= {x eD : dxix) <0}n{x e D : 0{x) ^ u} 

= {x e D -.iix) {0,u)} 
^ ^ (0,u)). 

Hence, for every q Cz denoting by j^q"'"'''^^'^^ the homology homomorphism in- 
duced by the inclusion D{$ ^ (a, u)) D{S ^ with (a,w) ^ (/?, w), it 
holds that 

PiK,0,^)Au,v) = rank (im7r(0'")'("^^^) . 
We claim that there exists a positive real number /3 such that 

for every /3 with < /3 < /3 (the claim is trivial for /3 = 0). In particular, this fact 
proves that P{K,<pfK)A'^^ ^) lim/3^o+ P{D,$),q ((0' (/^' ^)) ■ 
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In order to prove this claim, we consider the inverse system of homomorphisms 
^(o,u),ip,v) . ^ (o,w)) ^ Hq{D{^ < {13, V)) over the directed set {/3 e R : 

/3 > 0} decreasingly ordered. The foUowing isomorphisms hold: 

ij^^(o.S).(o,?) ^ -^^ - hmim^(°'")'(''^^T. 

Indeed, imTTg*^' = im ^imTTg"'"''^^'"'' by the continuity of Cech homology, and 

im limTTg^'"^'^^'^'* = lim imTr^'^'"^'*'^'^^ because the inverse limit of vector spaces is 
an exact functor and therefore it preserves epimorphisms, and hence images. 

It remains to prove that there exists a positive real number /3 such that, for 
every < /3 < /3, imTTg'^'"-''*-'^'''^ is isomorphic to ^unimTTg^'"^''-'^'^-'. To this end, let 
us consider the following commutative diagram, with < /?' < /3": 



H,{D{^<{Q,u))) ^H,{D{<^^{Q,u))) 



(2.1) 



^(o,a),(f)',B) 



^(0,it),(/3",i?) 



^(P' ,v),(li" ,v) 



From the above diagram (|2.ip . we see that each Wq^ ^^)'(^ induces a map Tq^ 

surjective. On the other hand, by the finiteness of the rank of im tt^'^'"'''^*''^'' and 
the monotonicity of the rank invariants, there exists /? > such that the rank of 
^(o,«),(/3 ,v) g]-^j^g Q^ij^ equal to the rank of TTq^'"^'^'' whenever < /?' < /?" < /3. 
Hence the maps rjj^ ^ are surjections between vector spaces of the same finite 
dimension, i.e. isomorphisms for every < /3' < /3" < ^. Thus, ^imimTrg'^'"''^^'^^ 
is the inverse limit of a system of finite dimensional vector spaces isomorphic to 

(0,ir),(^,'!7) • ,1 1 ■ 

im TTq , proving tlie claim. 

We now claim that for every strictly positive real number /3, there exists a strictly 
positive real number a < f3 such that 

for every a with < a < d. 

This claim can be proved in much the same way as the previous one. We consider 
the inverse system of homomorphisms TTg"'"'''^'"'' : Hq{D{$ ^ ia,u)) Hg{D{$ ^ 
{P,v}) over the directed set {a € M : < a < /?} decreasingly ordered. The 
following isomorphisms follow again from the continuity of Cech homology and the 
exacteness of the inverse limit functor for vector spaces: 

im7r(°'")'('''^") ^ im lim tt^"'")'''''^) limimTr^"'")'^'''"). 



imTTq"'"''^^ —J- imTTq"'"'''^'^ '"\ From diagram (|2.ip we see that these maps are 



To prove that there exists a strictly positive real number a such that, for every 

im 

following commutative diagram, with < a' < a" : 



< q; < d, imTTq"'"'*'*-''''''' is isomorphic to hmimTTg"'"''''''^''''', let us consider the 
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(2.2) (c',a),(/3.if) 



H,iD{<i> < {a',u))) ^H,{D{<f ^ {a",u))) 



H,{D{^ ^ v))) H,iD{^ ^ (/3, v))). 



From the above diagram (|2.2p . we see that each tt^" induces a map cr^" 

imTTg" imTTg" ^")'('''^)^ From diagram (|2.2p we see that these maps are 

injective. On the other hand, by the finiteness of the rank of imTTg"'"^'''^'^'', for any 
a with < a < /?, and the monotonicity of the rank invariants, there exists d, 
with < d < /3, such that the rank of tt," jg finite and equal to the rank 

of TTg" '")^(^'")^ whenever < a' < a" < a. Hence the maps aq" are injections 
between vector spaces of the same finite dimension, i.e. isomorphisms for every 
< a' < a" < a. Thus, limimTTg"'"'''^^'"^ is the inverse hmit of a system of finite 

dimensional vector spaces isomorphic to im 7rg"'"^'^'''"\ proving the claim. □ 

Many applications require that the presence of single outliers does not affect 
the evaluation of similarity. In these cases, always assuming K triangulable, it 
is sufficient to study the closure of the interior of K instead of K itself. Indeed, 
applying Theorems 12.11 and 12.21 with the closure of the interior of K instead of 
K, we obtain a result of stability of persistent homology groups with respect to 
the perturbations of the studied set and a reconstruction result for the original 
persistent homology groups modulo perturbations of zero measure. 

We underline once more that the results of this section are based on the idea of 
translating the problem of stability with respect to set perturbations into that of 
stability with respect to function perturbations. Therefore, the use of the distance 
function is only one among many ways to achieve this end and has the advantage 
of working well when sets are compared using the Hausdorff distance. One could 
conceive different ways, in connection with other methods to compare sets, as the 
following sections show. 

3. Stability with respect to other distances between sets 

Our approach can be easily adapted to noise that is small with respect to dis- 
tances other than the Hausdorff distance 6h- 

We first show how persistent homology can be made stable with respect to per- 
turbations of the sets measured using the symmetric difference distance (Theorem 
13. Then we show the stability with respect to perturbations of fuzzy sets (The- 
orem 13.21) . 



3.1. Stability with respect to the symmetric difference distance. We work 
with a non-empty closed subset K of a triangulable set D in R". In this case, 
instead of the distance function (Ik, our construction depends on the use of functions 
X% : ^- R, with e € M, e > 0, defined as 



X^Kix) ^ ^l{B,)-' ■ / XK{y)dy 
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where B^(x) denotes the n-disk centered at x with radius e, ~ 5^(0), and xk 
denotes the characteristic function of K. The underlying idea of this choice is that 
the closer a point x of the real plane is to a large part of the closer the value 
of (a;) is to 1. More precisely, \%;{x) = 1 if and only if ijl{B^{x) n K) — fi{B^), 
whereas Af^ (x) = if and only if fi{B^{x) n K) = 0. Clearly, Af^ is a continuous 
function for every real number e > 0. 

Analogously to Theorem 12.11 in this case we have the following result. 

Theorem 3.1. Let Ki,K2 be non-empty closed subsets of a triangulable subspace 
D o/M". Moreover, let 0i,02 '■ D M.^ be vector-valued continuous functions. 
Then, defining ^f,*| : D E'^+i by *f = (-Aff^,(^i) and *| = (-A^^,(^2), the 
following inequality holds: 

(3.1) D^atch(p(D,^i),q^P{D,^i),^ <max -(^2j|«,|. 



Proof. For every x £ D, 



XkAv) - XK-Av) dy 

yeB^(x) 

XkAv) - XK^iy)] dy 



= ll{Be)-^ ■ pi{KiAK2). 

Thus \\\%^^ - AffJIoo < KBeY^ ■ n{KiAK2). The Muhidimensional Stability The- 
orem 11.51 for measuring function perturbations implies that 

Dmatch - 11*1 ^ *2l|oo- 

It follows that 

Drnatch(^P(^D,^l),q^P{D,^f,),q) ^ "lax { 1 1 A|f ^ - A|-J j oo , 1 1 " <^2 1 1 oo } 

< max{^(Be)"^ • fi{KiAK2), \\(pi - <^2||oo} 
= max{^(Be)"^ • dA{Ki,K2), - 02\\oo] ■ 



□ 



The previous theorem shows that, under our hypotheses, if two compact subsets 
Ki , K2 of the real plane are close to each other in the sense that their symmetric 
difference has a small measure, then also the rank invariants constructed using the 
functions , ^1 ^I'S close to each other. 

We observe that the estimate in inequality p.ip can be improved by substituting 
d-A {Ki , K2) with max^gR,> /^^^^ 



XKi(y) - XK2{y) dy 



3.2. Stability with respect to perturbations of fuzzy sets. Now we consider 
the case when sets are defined according to fuzzy theory, that is through functions 
representing the grade of membership of points to the considered set. One obtains 
a fuzzy set, for example, when a probability density p{x) is given, p{x) express- 
ing the probability that a point of the considered set belongs to an infinitesimal 
neighborhood of x. We confine ourselves to considering only probability densities 
with compact support contained in a triangulable subspace D of M". From the 
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Multidimensional Stability Theorem 11.51 for measuring function perturbations we 
immediately deduce the following result, whose simple proof is omitted, concern- 
ing the stability with respect to perturbations of fuzzy sets defined by probability 
densities. 

Theorem 3.2. Letpi,p2 be two probability density functions having support con- 
tained in a compact and triangulable subspace D of W . Defining ^'i,5'2 : D — > 
]gfc+i vj/j^ — (— and ^2 = {^P2t02), the following statement holds: 

4. An example 

In this section the theoretical framework presented in Section [5] is applied in 
a discrete setting. Our goal is to check the stability of the proposed framework 
with respect to set perturbations measured by the Hausdorff distance. We confine 
ourselves to the case q = 0. 

With this aim in mind, we work with the binary digital image represented in 
Figure [1] (left), and we corrupt this image by adding salt & pepper noise to a 
neighborhood of the set of its black pixels, as shown in Figure [1] (right). 




Figure 1. Two binary images of an octopus. The image on the 
right is a noisy version of that on the left. 

Black pixels of left and right images represent the sets Ki, K2 under study, 
respectively, whereas in both cases the 269x256 rectangle of black and white pixels 
together constitute the set D. The so obtained noisy set K2 is close to the original 
set Ki with respect to the Hausdorff distance. 

A graph structure based on the local 4-neighbors adjacency relations of the 
digital points is used in order to topologize the images. 

Fixed the point c G D corresponding to the center of mass of Ki, the chosen 
measuring function for both instances is 1^ : _D M, (p{p) = —\\p — c\\. 

Figure [2] (left) shows the persistence diagram of the 1-dimensional 0th rank 
invariant P(Ki,ip^Ki},o- displays eight relevant points in the persistence diagram, 
corresponding to the eight tentacles of the octopus. Only one of these points is 
at infinity (and therefore depicted by a vertical line rather than by a circle) since 
Ki has only one connected component. As for P(^K2,ip\K2),o^ presence 
of a great quantity of connected components in the noisy octopus, its persistence 
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diagram has a very large number of points at infinity, and a figure showing them 
all would be hardly readable. For this reason Figure [2] (right) shows only a small 
subset of its persistence diagram. However it is sufficient to perceive how dissimilar 
it is from P(ifi,^|^j,o- 




Figure 2. Left: The persistence diagram of the rank invariant 
P{Ki,ip^K-^},o corresponding to the original octopus image. Right: A 
detail of the persistence diagram of the rank invariant P(k,,i^|j^^),o 
corresponding to the noisy octopus image. 



As suggested by Theorem 12. li if instead we compare Ki and K2 by means of 
the rank invariants ^ and p^^ ^^-j q, where $1 : D — R^, $1 = (d-Ki , and 

$2 : — K^, $2 = {dK2i </')i can see the similarity between Ki and K2 modulo 
the salt & pepper noise. This is illustrated in Figure[3l In Figure[3](a)-(b), we show 
the rank invariants p^j^ q and p^jj ^ both restricted to the half-plane '''(fg), 

with r= (0.1483,0.9889) and 6 = (13.0434,-13.0434), that is the half-plane of the 
foliation containing the point ((0, —100), (3, —80)). In other words. Figure [3] (a)- (b) 
shows p J J and p 4^ , respectively. We can appreciate their similarity, 

even if their matching distance dmatch is not necessarily smaller than the Hausdorff 
distance between Ki and i^2- The considered half-plane has been chosen so that it 
contains points where the rank invariant takes non-trivial values. 

Indeed it is easy to verify that Theorem 12.11 does not guarantee the stability of 
P(D F* 1 stability of , , where p = miu; li. We point out that 

pf { N {p ■ p ■ t) = p/ 4 \ (s, t), and hence the passage from the mea- 

suring function -F,*-, to the measuring function p ■ corresponds to "rescaling 
up" the domain of the rank invariant. In other words, when we change Ki into a 
new compact set K2 that is close to Ki with respect to the Hausdorff distance, the 
matching distance between P{Ki,ip\Ki),q ^^"^ P{K-,,ip\K2),q '^'^^ small, while the 

one between 0/ (d„ ^.^ )\ and 0/ (d^- ^.^ )\ must be small. 
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This is illustrated in FiffureEl where the rank invariants «, . andp, , 

^ {l.h)'' ^ ' (I.b)-"' 

displayed the top row, are not as similar as the rank invariants , and 

p it,^ , displayed in the bottom row. 

^ (1,6)'' 




-95 -90 -85 -80 -75 -70 -65 -60 -55 -50 



(a) 




-20 -15 -10 -5 5 10 15 20 25 



fl ■ s 




-95 -90 -85 -80 -75 -70 -65 -60 -55 -50 



(b) 




-20 -15 -10 -5 5 10 15 20 25 



fl ■ s 



(c) (d) 

Figure 3. (a) The rank invariant p^j^ ^ restricted to the half- 
plane TT^j-^y with r= (0.1483,0.9889) and = (13.0434,-13.0434), 
that is the half-plane of the foliation containing the point 
((0, —100), (3, —80)). (b) The rank invariant P(£) q restricted 
to the same half-plane, (c)-(d) The same restrictions as in (a)-(b), 
respectively, but rescaled by ^ = min{?i,/2}- 



Next we show how it is possible to point-wisely recover the rank invariant of 
{Ki,ip\Xi) from that of (£>, $i). According to Theorem 12.21 P(_r-i.^|^,j.o(^': w) = 
P{D <ii)('^' 1^' ^) ^'-'^ a, /? > sufficiently small. 
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As shown in [I], in this case the fohowing equahties hold (with reference to 
Definition [Q]): 

h{a,u,l3,v)= =, l2{a,u,p,v)^ =, 

(4.1) bAa u 3 v) - u 3 v\ - 

s(a,u,/3,z;) = i^, /j, ^) = 

As a consequence, Theorems 12.21 and 11.41 (appUed in this order) imply that for 
every pair (u, v), with u < v, and for < a < (3, with a and 13 sufficiently small, 



Y 



p/ 5 X (s(a,u,^,?;),i(a,'u,/3,w)) 

D.F J-^ .q 



a + li /3 - 



where : D ^ M is defined by setting, for every x € D, 



Hence the finite value P(Ki,y|^_i),g(w, u) is equal to p, x ^ (^ffpf;, 17^), 



if we choose a and /3 small enough. The corresponding admissible pair (Z, 6) results 
to be close to the pair ((0, 1), (0, 0)). 

In other words, the information about the rank invariant of the original pair 
(-ftTi, '^\K-^ can be recovered on the leaves associated with the admissible pairs (/, 6) 
in a small neighborhood of the pair ((0, 1), (0, 0)), after re-paramcterizing these 

leaves by taking each point (m, u) to the point if+jj)- Note that the pair 

((0, 1), (0,0)) is not admissible but is located on the boundary of the set Adm-i. 
This leads to instabilities if wc take a, /3 too small. 

We underline that ThGorem l2.2l cnsures this approximation is good only point-by- 
point. Thus, even for admissible pairs (I, b) close enough to the pair ((0, 1), (0, 0)), 
the matching distance between p/ ^ \ and P(K,ip^K },q quite large. 

To illustrate these issues we have computed the value taken by P(^Ki,'p\k^),o at the 
point (u,w) = (—100,-80), that is 8. Using formulas (|4.ip . we have computed the 
admissible pairs / = /(a,/3), b = &(a, /3) such that the half-plane contains the 
point ((a, —100), (/3, —80)) for the values of a,/3 shown in Table 111 For the same 
values of a,/?, Table |3] also shows the parameters s = s(a, /3) and t — t{a, (3) for 
which we have that P(^^^j q((q;, -100), (/3, -80)) = p/ ^ \ (s,t). 

\^'^(r.b) J 

Computations show that J.o(^100, — 80) = Pf \ (s,i) = 8 for 



small but positive values of a and /3. It is noticeable that for a = 0, due to the 
mentioned instabilities near the boundary of the set Adm2, computations are not 
reliable. 
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U 


p 


V 


s 


+ 

L 


n , s ( K i\ 


0.5 


-100 


24 


-80 


-70.5866 


-39.7272 


1 


0.5 


-100 


16 


-80 


-70.9216 


-45.6179 


3 


0.5 


-100 


8 


-80 


-77.2843 


-55.9262 


3 


0.5 


-100 


1 


-80 


-97.1120 


-77.1040 


8 


0.5 


-100 


0.65 


-80 


-98.7692 


-78. .7672 


8 


0.3 


-100 


0.45 


-80 


-98.9677 


-78.9657 


8 


0.1 


-100 


0.25 


-80 


-99.1663 


-79.1643 


8 





-100 


0.15 


-80 


-99.2655 


-79.2635 






Table 1. The parameters used to approximate the value of 
P{KuviK^},o at {u,v) = (-100,-80), that is 8, using P f \ ■ 



The corresponding rank invariants p / \ for the choices of a and /3 

considered in Table S] are displayed in Figure SI 



5. Discussion 

In this paper we have shown the stability of persistent homology groups with 
respect to perturbations of the studied set. Measuring set perturbations by different 
distances requires different constructions in order to achieve stability. 

If set perturbations are measured through the Hausdorff distance, we replace 
compact sets by distance functions. In this way, by the well-known property that if 
K' is a good Hausdorff approximation of K then the distance function dx' is close 
to dx, and by utilizing already available results of persistent homology stability 
with respect to function perturbations, we deduce stability with respect to set 
perturbations. We also show that while passing from sets to functions we are still 
able to recover information about the persistent homology groups of the original 
sets. 

If set perturbations arc measured through the symmetric difference distance, an 
analogous procedure leads to the proof of stability also in this case. Finally, using 
the sup-distance enables us to guarantee stability also with respect to perturbations 
of fuzzy sets. 

The common underlying idea is to compare sets by comparing functions describ- 
ing the sets themselves. 

While considering the Hausdorff distance and the symmetric difference distance 
is certainly not exhaustive of all the possible ways of measuring set perturbations, 
it accounts for two very widely used ones, and allows us to indicate a general 
procedure that could be applied also when dealing with other distances. 

Finally, we underline that the technique developed in this paper essentially relies 
on the multidimensional generalization of persistent homology, showing once more 
the importance of further pursuing this area of research. 
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Figure 4. The rank invariant p / , \ as a and /3 tend 

to 0. Red circles and red lines denote the points (proper or at infin- 
ity) of the corresponding persistence diagram. The blue diamonds 
denote the point (s,t) corresponding to ((a, —100), —80)). 



